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Polarization measurements are the best instrument to understand how quark and an- 
1 1 tiquark combine into the different quarkonium states, but no model has so far succeeded 

in explaining the measured J/tp and T polarizations. On the other hand, the experimen- 
tal data in proton- ant iproton and proton-nucleus collisions are inconsistent, incomplete 
and ambiguous. New analyses will have to properly address often underestimated issues: 
the existence of azimuthal anisotropics, the dependence on the reference frame, the in- 
fluence of the experimental acceptance on the comparison with other measurements and 
with theory. Additionally, a recently developed frame-invariant formalism will provide 
an alternative and often more immediate physical viewpoint and, at the same time, will 
help probing systematic effects due to experimental biases. The role of feed-down de- 
cays from heavier states, a crucial missing piece in the current experimental knowledge, 
^v^j will have to be investigated. Ultimately, quarkonium polarization measurements will also 

I offer new possibilities in the study of the properties of the quark-gluon plasma. 
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1. The experimental situation 

Quarkonia, bound states made of a quark (of type charm or beauty) and its anti- 
quark, offer us a privileged window over the physics of the strong force, which is at 
the origin of visible matter and, yet, is the least well-understood aspect of the Stan- 
dard Model of elementary interactions. Quarkonia represent the most elementary 
manifestation of the strong binding force and allow us to study crucial open ques- 
tions: how are quarks confined inside hadrons? How do strong forces generate the 
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properties of particles made of quarks? Can quarks become unbound under extreme 
conditions (high temperature and density: the quark-gluon plasma), as they existed 
in the first moments of the universe? To test and consolidate the current theory 
of the strong force, quantum chromo-dynamics (QCD), it is crucial to study how 
quarkonia are produced in elementary (proton-proton) collisions and in the much 
more complex nucleus-nucleus collisions, where the potential that binds the quarks 
and the gluons should be screened and the medium should reflect the partonic de- 
grees of freedom. However, our present understanding of this physics topic is rather 
limited, despite the multitude of experimental data accumulated over more than 30 
yearsPThe J/ip and ip f direct production cross sections measured (in the mid 1990's) 
by CDF, in pp collisions at 1.8 TeV, 2 were seen to be around 50 times larger than 
the available expectations, based on leading order calculations made in the scope of 
the Colour Singlet Model (CSM). The non-relativistic QCD (NRQCD) framework, 3 
where quarkonia can also be produced as coloured quark pairs, succeeded in describ- 
ing the measurements, opening a new chapter in the studies of quarkonium produc- 
tion physics. However, these calculations depend on non-perturbative parameters, 
the long distance colour octet matrix elements, which have been freely adjusted to 
the data, thereby decreasing the impact of the resulting agreement between data 
and calculations. More recently, calculations of next-to-leading-order (NLO) QCD 
corrections to colour-singlet quarkonium production showed an important increase 
of the high-px rate, significantly decreasing the colour-octet component needed to 
reproduce the quarkonium production cross sections measured at the Tevatron. 4 
Given this situation, differential cross sections are clearly insufficient information 
to ensure further progress in our understanding of quarkonium production. Polar- 
ization measurements, determining the average angular momentum states of the 
produced quarkonia from their decay distributions, can provide the definitive tests 
of the theory of quarkonium production. No other study addresses more directly the 
question: how does the observed quark-antiquark bound states acquire their final 
quantum numbers? In fact, the competing mechanisms dominating in the different 
theoretical approaches lead to very different expected polarizations of the produced 
quarkonia at high p T . NRQCD predictdMQ almost fully transverse polarization 
(angular momentum component J z = ±1) for directly produced ^ and J/i/j mesons 
with respect to their own momentum direction (the helicity frame), while according 
to the new NLO calculations of colour-singlet quarkonium production^ these states 
should show a strong longitudinal (J z = 0) polarization component. 

Having two very different theoretical predictions appears to be an ideal situ- 
ation in the prospect of discriminating between the two theory frameworks using 
experimental data. However, the present experimental knowledge is incomplete and 
contradictory. A significant fraction (around one thircP) of promptly produced J/ip 
mesons (i.e. excluding contributions from B hadron decays) comes from \c and ip' 
feed-down decays. This sizeable source of indirectly produced J/ifj mesons is not 
subtracted from the current measurements, and its kinematic dependence is not 
precisely known. Despite this limitation, it seems safe to say that the pattern mea- 
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sured by CDF^of a slightly longitudinal polarization of the inclusive prompt J ftp is 
incompatible with any of the two theory approaches mentioned above. The situation 
is further complicated by the intriguing lack of continuity between fixed-target and 
collider results, which can only be interpreted in the framework of some specific 
(and speculative) assumptions still to be tested. 10 

The bb system should satisfy the non-relativistic approximation much better 
than the cc case. For this reason, the T data are expected to represent the most 
decisive test of NRQCD. However, the present data from Tevatron) 11 ! 12 ! 13 ! for 
(Pt) < 27 GeV/c, tend to contradict the crucial NRQCD hypothesis that high-px 
quarkonia, produced by the fragmentation of an outgoing (almost on-shell) gluon, 
are fully transversely polarized along their own direction. At lower energy and px? 
the E866 experiment 14 has shown yet a different polarization pattern: the T(2S') and 
T(35') states have maximal transverse polarization, with no significant dependence 
on transverse or longitudinal momentum, with respect to the direction of motion of 
the colliding hadrons (Collins-Soper frame^). Unexpectedly, the Y(IS'), whose spin 
and angular momentum properties are identical to the ones of the heavier T states, 
is, instead, found to be only weakly polarized. These results give interesting physical 
indications. First, the maximal polarization of T(2S) and T(3S) along the direction 
of the interacting particles places strong constraints on the topology and spin prop- 
erties of the underlying elementary production process. Second, the small T(1S) 
polarization suggests that the bottomonium family may have a peculiar feed-down 
hierarchy, with a very significant fraction of the lower mass state being produced in- 
directly; at the same time, the polarization of the T's coming from Xb decays should 
be substantially different from the polarization of the directly produced ones. 

This rather confusing situation demands a significant improvement in the ac- 
curacy and detail of the polarization measurements, ideally distinguishing between 
the properties of directly and indirectly produced states. We remind that the lack 
of a consistent description of the polarization properties represents today's biggest 
uncertainty in the simulation of the LHC quarkonium production measurements 
and represents the largest contribution to the systematic error affecting the mea- 
surements of quarkonium production cross sections and kinematic distributions. 

It is true that measurements of the quarkonium decay angular distributions are 
challenging, multi-dimensional kinematic problems, requiring large event samples 
and a very high level of accuracy in the subtraction of spurious kinematic corre- 
lations induced by the detector acceptance. The complexity of the experimental 
problems which have to be faced in the polarization measurements is testified, for 
example, by the disagreement between the CDF J/ip results obtained in Run 
and Run II 9 and by the contradictory T(1S) results obtained by CDF^and DCP^. 
However, it is also true, as we shall emphasize hereafter, that most experiments 
have presented in the published reports only a fraction of the physical information 
derivable from the data. This happens, for example, when the measurement is per- 
formed in only one polarization frame and is limited to the polar projection of the 
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decay angular distribution. As we have already argued in Ref. [101 these incomplete 
measurements do not allow definite physical conclusions. At best, they confine such 
conclusions to a genuinely model-dependent framework. Moreover, such a fragmen- 
tary description of the observed physical process obviously reduces the chances of 
detecting possible biases induced by not fully controlled systematic effects. 

In this work we focus our attention on aspects that need to be taken in consider- 
ation in the analysis of the data, so as to maximize the physical significance of the 
measurement and provide all elements for its unambiguous interpretation within 
any theoretical framework (Sects. (2jj5|. 

We also discuss (Sect. [6| how the polarization of vector quarkonia, measured 
from dilepton event samples, can be used as an instrument to study the suppression 
of Xc and Xb in heavy-ion collisions, where a direct determination of signal yields 
involving the identification of low-energy photons is essentially impossible. 

2. Basic concepts 

Because of angular momentum conservation and basic symmetries of the electro- 
magnetic and strong interactions, a particle produced in a certain superposition 
of elementary mechanisms may be observed preferentially in a state belonging to 
a definite subset of the possible eigenstates of the angular momentum component 
J z along a characteristic quantization axis. When this happens, the particle is said 
to be polarized. In the dilepton decay of quarkonium, the geometrical shape of 
the angular distribution of the two decay products (emitted back-to-back in the 
quarkonium rest frame) reflects the average polarization of the quarkonium state. 
A spherically symmetric distribution would mean that the quarkonium would be, 
on average, unpolarized. Anisotropic distributions signal polarized production. 

The measurement of the distribution requires the choice of a coordinate system, 
with respect to which the momentum of one of the two decay products is expressed 
in spherical coordinates. In inclusive quarkonium measurements, the axes of the 
coordinate system are fixed with respect to the physical reference provided by the 
directions of the two colliding beams as seen from the quarkonium rest frame. The 
polar and azimuthal angles d and tp describe the direction of one of the two decay 
products (e.g. the positive lepton) with respect to the chosen polar axis and to the 
plane containing the momenta of the colliding beams ("production plane"). The 
actual definition of the decay reference frame with respect to the beam directions 
is not unique. Measurements of the quarkonium decay distributions used mainly 
two different conventions for the orientation of the polar axis: the flight direction 
of the quarkonium itself in the centre-of-mass of the colliding beams (centre-of- 
mass helicity frame, HX) and the bisector of the angle between one beam and the 
opposite of the other beam (Collins-Soper frame, CS). The motivation of the latter 
definition is that, in hadronic collisions, it coincides with the direction of the relative 
motion of the colliding partons, when their primordial transverse momenta, &t, are 
neglected. We note that these two frames differ by a rotation of 90° around the y axis 
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when the quarkonium is produced at high p T and negligible longitudinal momentum 
(Pt ^> |Pl|)- All definitions become coincident in the limit of zero quarkonium • In 
this limit, moreover, for symmetry reasons any azimuthal dependence of the decay 
distribution is physically forbidden. 

The most general decay angular distribution for inclusively observed quarkonium 
states can be written as^3 

W(cos oc - — r (1 + A# cos 2 $ + sin 2 $ cos 2ip + A#<^ sin 2$ cos cp) , (1) 

(3 + A#) 

where the three parameters A# ("polarization"), X^ and A#<^ satisfy the relations^ 
|A^| < |(l + A^), A 2 +2A 2 ^<1, 

|A*„|<^(1-A„), (2) 

(1 + 2A^) 2 + 2A^ < 1 for A^ < -1/3 , 

which, in particular, imply |A^| < 1, |A#^| < a/2/2, |A^| < 0.5 for A# = and 
A^ -> for A# — » -1. 

3. The importance of the reference frame and of the azimuthal 
anisotropy 

The coefficients A$ , A^ and A^^ depend on the polarization frame. To illustrate the 
importance of the choice of the polarization frame, we consider specific examples 
assuming, for simplicity, that the observation axis is perpendicular to the natural 
axis. This case is of physical relevance since when the decaying particle is produced 
with small longitudinal momentum (\ph\ <^ Pt, a frequent kinematic configuration 
in collider experiments) the CS and HX frames are actually perpendicular to one 
another. In this situation, a natural "transverse" polarization (A# = +1 and A^, = 
^$ip = 0), for example, transforms into an observed polarization of opposite sign 
(but not fully "longitudinal"), A^ = —1/3, with a significant azimuthal anisotropy, 
X'^ = 1/3. In terms of angular momentum wave functions, a state which is fully 
"transverse" with respect to one quantization axis (|J, J z ) = |1,±1)) is a coherent 
superposition of 50% "transverse" and 50% "longitudinal" components with respect 
to an axis rotated by 90°: 

|1,±1> ^ \ + \ T ^= |1,0). (3) 

The decay distribution of such a "mixed" state is azimuthally anisotropic. The same 
polar anisotropy A^ = —1/3 would be measured in the presence of a mixture of at 
least two different production processes resulting in 50% "transverse" (|J, J z ) = 
1 1 , =bl) ) and 50% "longitudinal" (\J,J Z ) — 1 1 , 0) ) natural polarization along the 
chosen axis. In this case, however, no azimuthal anisotropy would be observed. As 
a second example, we note that a fully "longitudinal" natural polarization (A# = 
— 1) translates, in a frame rotated by 90° with respect to the natural one, into a 



August 24, 2012 1:0 WSPC/INSTRUCTION FILE QuarkoniumPolariza- 
tion 



6 Pietro Faccioli 

fully "transverse" polarization = +1), accompanied by a maximal azimuthal 
anisotropy (A^ = —1). In terms of angular momentum, the measurement in the 
rotated frame is performed on a coherent admixture of states, 

11,0} ^ -J=|l,+1>-^=|1,-1), (4) 

while a natural "transverse" polarization would originate from the statistical su- 
perposition of uncorrelated |1,+1) and |1,— 1) states. The two physically very dif- 
ferent cases of a natural transverse polarization observed in the natural frame and 
a natural longitudinal polarization observed in a rotated frame are experimentally 
indistinguishable when the azimuthal anisotropy parameter is integrated out. These 
examples show that a measurement (or theoretical calculation) consisting only in 
the determination of the polar parameter A# in one frame contains an ambiguity 
which prevents fundamental (model-independent) interpretations of the results. The 
polarization is only fully determined when both the polar and the azimuthal com- 
ponents of the decay distribution are known, or when the distribution is analyzed 
in at least two geometrically complementary frames. 

Due to their frame-dependence, the parameters A#, and A#<^ can be affected 
by a strong explicit kinematic dependence, reflecting the change in direction of 
the chosen experimental axis (with respect to the "natural axis") as a function 
of the quarkonium momentum. As an example, we show in Fig. [I] how a natural 
transverse J/i/j polarization (A# = +1) in the CS frame (with A^ = A$^ = and no 
intrinsic kinematic dependence) translates into different px-dependent polarizations 
measured in the HX frame in different rapidity acceptance windows, representative 
of the acceptance ranges of several Tevatron and LHC experiments. This example 




P T [GeV/c] p T [GeV/c] 



Fig. 1. Kinematic dependence of the J/ip decay angular distribution seen in the HX frame, for a 
natural polarization \$ = +1 in the CS frame. The curves correspond to different rapidity intervals; 
from the solid line: \y\ < 0.6 (CDF), \y\ < 0.9 (ALICE), \y\ < 1.8 (DO), \y\ < 2.5 (ATLAS and 
CMS), 2 < y < 5 (LHCb). For simplicity the event populations were generated flat in rapidity. 

shows that an "unlucky" choice of the observation frame may lead to a rather 
misleading representation of the experimental result. Moreover, the strong kinematic 
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dependence induced by such a choice may mimic and/or mask the fundamental 
("intrinsic") dependencies reflecting the production mechanisms. 
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Fig. 2. Polar anisotropy of the J/ip decay distribution as seen in the CS (left) and HX (right) 
frames, when all the events have full transverse polarization, but 60% in the CS frame and 40% 
in the HX frame. The curves represent measurements in different rapidity ranges (see Fig. [T|. 

Not always an "optimal" quantization axis exists. This is shown in Fig. [2] where 
we consider, for illustration, that 60% of the J/i/j events have natural polarization 
A# = +1 in the CS frame while the remaining fraction has A# = +1 in the HX frame. 
Although the polarizations of the two event subsamples are intrinsically indepen- 
dent of the production kinematics, in neither frame, CS or HX, will measurements 
performed in different transverse and longitudinal momenta windows find "simple" , 
identical results. Corresponding figures for the T(1S) case can be seen in Ref. fTUI 
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Fig. 3. The CDF J/ip polarization measurement in the helicity frame (data points) and the range 
for the corresponding polarization in the CS frame (dashed line: CS polarization for A^ x = 0). 



CDF measured for the J/ip an almost vanishing polar anisotropy parameter in 
the helicity frame. It is natural to wonder how the measurement would look like 
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in the CS frame. However, the transformation to another frame depends on the 
azimuthal anisotropy, which was not reported by the experiment. For example, as 
shown in Fig. [3j if the distribution in the HX frame were azimuthally isotropic, 
the measured polarization would correspond to a practically undetectable polariza- 
tion in the CS frame (dashed line). However, if we take into account all physically 
possible values of the azimuthal anisotropy, as allowed by the relations in Eq. [2| 
we can only derive a broad spectrum of possible CS polarizations, approximately 
included between —0.5 and +1 (shaded band). This example shows how a measure- 
ment reporting only the polar anisotropy is amenable to several interpretations in 
fundamental terms, often corresponding to drastically different physical cases. 

An analysis ignoring the azimuthal dimension can produce wrong results. In 
fact, the experimental acceptances for the variables cos$ and cp are usually strongly 
intercorrelated because of the limited sensitivity to low- momentum leptons, which 
reduces the population of events in specific angular regions, depending on the ref- 
erence frame. For example, the experimental efficiency for the projected cos^ dis- 
tribution depends on the ip distribution, that is on A^, and vice- versa. If the cp 
dimension is integrated out and ignored, the A# measurement becomes strongly 
dependent on the specific "prior hypothesis" (implicitly) made for the angular dis- 
tribution in the Monte Carlo simulation. To illustrate this concept, we consider J/ip 
pseudo-data in the kinematic region 9 < pt < 12 GeV/c, < \y\ < 1, simulating 
the acceptance filter with the requirement that both leptons have pt > 3 GeV/c. 
The angular acceptances for these conditions in the CS and HX frames are shown 
in Fig. [4j We consider the example scenario of a fully longitudinal polarization 
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Fig. 4. Angular acceptances in the CS and HX frames for J/ip decay into leptons, in the kinematic 
region 9 < px < 12 GeV/c, < \y\ < 1, when only leptons having > 3 GeV/c are detected. 



in the HX frame. A one-dimensional measurement is performed in the CS frame 
integrating out and ignoring the cp dependence. The detector- acceptance correc- 
tion is performed one-dimensionally, using Monte Carlo data generated assuming 
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Fig. 5. Results of a pseudo-experiment (~ 60k reconstructed dilepton events) where the J/ip po- 
larization (generated as fully longitudinal in the HX frame) is measured through one-dimensional 
angular projections in the CS and HX frames. a,b,e,f: (wrong) results obtained using a "standard" 
unpolarized Monte Carlo simulation for the acceptance correction. c,d,g,h: ("correct") results ob- 
tained after reweighing iteratively the Monte Carlo data according to the results found. 
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a flat azimuthal dependence. Figure [5^i shows that the acceptance-corrected cost? 
distribution in the CS frame is flat, leading to a wrong "unpolarized" result, re- 
flecting the polarization assumption used in the Monte Carlo simulation. If the 
Monte Carlo data used for the acceptance correction are reweighted to the "true" 
polarization (a two-dimensional ingredient), the same distribution changes drasti- 
cally (Fig. [5^), correctly showing a strong transverse polarization (the CS frame 
being almost perpendicular to the HX frame). This shows that when only a one- 
dimensional projected distribution is measured, the detector acceptance description 
must, nevertheless, be maintained multi-dimensional. One-dimensional acceptance 
corrections or "template" fits should be avoided, unless the MC is iteratively re- 
generated with the correct distribution of the variables that have been integrated 
out (which has, therefore, to be measured anyway). Unfortunately, one-dimensional 
polarization analyses are widespread, even, paradoxically, in precision tests of the 
Standard Model and searches for new physics. To measure the polar anisotropy of W 
decays or Drell-Yan production using template distributions (to account for accep- 
tance and efficiency) integrated over the azimuthal angle, as done in recent analyses 
reported by LHC experiments, may strongly bias the measurement towards the dis- 
tribution used to produce the Monte Carlo simulation. One analysis even imposes 
the absence of azimuthal anisotropics, assuming that the data are exactly described 
by the naive Born- level Drell-Yan angular distribution valid at = 0. New physics 
effects changing drastically the azimuthal anisotropy with respect to the expected 
one (assumed in the Monte Carlo) may be missed by this kind of analyses. 

4. A frame- invariant approach 

It can be shown that the combination of coefficients 

A = (A# + 3A^)/(1 - A^) (5) 

is independent of the polarization frame, The fundamental meaning of the frame- 
invariance of this quantity is discussed in Ref. [20j The determination of A is immune 
to "extrinsic" kinematic dependencies induced by the observation perspective and 
is, therefore, less acceptance-dependent than the standard anisotropy parameters 
A# , A^ and A#^. Referring to the example shown in Fig.[2j any arbitrary choice of the 
experimental observation frame will always yield the value A = +1, independently 
of kinematics. This particular case, where all contributing processes are transversely 
polarized, is formally equivalent to the Lam- Tung relation.^ The existence of frame- 
invariant parameters also provides a useful instrument for experimental analyses. 
Checking, for example, that the same value of an invariant quantity is obtained, 
within systematic uncertainties, in two distinct polarization frames is a non-trivial 
verification of the absence of unaccounted systematic effects. In fact, detector geom- 
etry and/or data selection constraints strongly polarize the reconstructed dilepton 
events, as shown in Fig. [4j Background processes also affect the measured polariza- 
tion, if not well subtracted. The spurious anisotropics induced by detector effects 
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and background do not obey the frame transformation rules characteristic of a 
physical J = 1 state. If not well corrected and subtracted, these effects will distort 
the shape of the measured decay distribution differently in different polarization 
frames. In particular, they will violate the frame-independent relations between the 
angular parameters. Any two physical polarization axes (defined in the rest frame 
of the meson and belonging to the production plane) may be chosen to perform 
these "sanity tests" . The HX and CS frames are ideal choices at high and mid 
rapidity, where they tend to be orthogonal to each other. At forward rapidity and 
low pt? the significance of the test can be maximized by using the CS axis and 
the "perpendicular helicity axis" which coincides with the helicity axis at zero 
rapidity and remains orthogonal to the CS axis at nonzero rapidity. Given that A is 
"homogeneous" to the anisotropy parameters, the difference A^ B ^ —M A ^ between the 
results obtained in two frames provides a direct evaluation of the level of systematic 
uncertainties not accounted in the analysis. 

To illustrate the application of the frame-independent formalism as a tool to 
spot problems in experimental data analyses, we refer again to the above-described 
J/ip pseudo-experiments. The result of the measurement performed with one- 
dimensional acceptance correction assuming unpolarized production is shown in 
Fig. [5^i,b,e,f, including, this time, both the polar and azimuthal projections, in the 
CS and in the HX frame. From the comparison of these four one-dimensional re- 
sults we derive, using Eq.|5j A cs ~ -0.32 and A HX ~ -0.93. The difference between 
the two values is an unequivocal signal of a mistake in the analysis. In fact, after 
reweighting the Monte Carlo with the "correct" polarization (Fig. [5]:,d,g,h), which 
can be iteratively inferred choosing as "generation frame" the one showing at each 
step the strongest polarization modulations (the HX frame in this case), both A 
values approach —1, as expected in this exercise. 

Incidentally, we note that the "true" A value is generally not included between 
the values found in two different reference frames. The best value of A in the pres- 
ence of not completely corrected systematic effects is not the average among dif- 
ferent frames and is best approximated by the value found in the frame showing 
the smallest acceptance correlations between cos$ and ip (HX frame in the above 
example, see Fig. [4|. It is, therefore, a priori not fully justified to impose the con- 
straint A cs = A HX in a fit of the angular distributions performed simultaneously in 
two frames, as done in a recent LHC analysis of J ftp polarization. 

Another example of utility of the invariant polarization parameter can be seen in 
Fig. [6j showing J/iJj polarization "measurements" in the CS and HX frames versus 
px- While the A# values seem to change significantly from one frame to the other, 
the two A^ patterns are very similar. This observation alerts to an experimental 
artifact in the data analysis. We can evaluate the significance of the contradiction 
by calculating the frame- invariant A variable in each of the two frames. For the case 
illustrated in Fig. [ij averaging the four represented pt bins, we see that A in the 
HX frame is larger than in the CS frame by 0.5 (a rather large value, considering 
that the decay parameters are bound between —1 and +1). In other words, the de- 
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Fig. 6. Example of data where the J/ip polarization measurements in the CS and HX frames 
(empty and filled symbols, respectively) are inconsistent with each other. 

termination of the decay parameters must be biased by systematic errors of roughly 
this magnitude. Given the puzzles and contradictions existing in the published ex- 
perimental results, as recalled in Section [I] the use of a frame- invariant approach to 
perform self-consistency checks, which can expose unaccounted systematic effects 
due to detector limitations and analysis biases, constitutes a non-trivial comple- 
mentary aspect of the methodologies for quarkonium polarization measurements. 

5. The role of the feed-down decays 

Many of the prompt J/t/j and T mesons produced in hadronic collisions result from 
the decay of heavier S- or P-wave quarkonia. However, the existing polarization 
measurements at collider energies make no distinction between directly and indi- 
rectly produced states. The role of the feed-down from heavier S states (responsible, 
for example, for about 8% of J/ip production at low p^ 8 ) is rather well understood. 
Data of the BES 23 and CLE0EH3 experiments in e + e collisions indicate that in 
the decays %p' — >> J/ifjirir and T(2S) —> T(1S)7T7t the di-pion system is produced 
predominantly in the spatially isotropic (5- wave) configuration, meaning that no 
angular momentum is transferred to it. Consequently, the angular momentum align- 
ment is preserved in the transition from the 2S to the IS state. This allows us to 
assume that the dilepton decay angular distribution of the J/ip [T(15)] mesons re- 
sulting from %p' [T(2/35)] decays is the same as the one of the %p' [T(2/35)], provided 
that a common polarization axis is chosen for the two particles. At high momentum, 
when the J/ip and ip' directions with respect to the centre of mass of the colliding 
hadrons practically coincide, ip' mesons and J/iJj mesons from ijj' decays have the 
same observable polarization with respect to any system of axes defined on the basis 
of the directions of the colliding hadrons. In the case of the polar anisotropy param- 
eter A#, for instance, the relative error, |AA#/A#|, induced by the approximation of 
considering the J/iJj and ip' directions as coinciding is (9[(Am/p) 2 ], where Am is 
the 2S — IS mass difference and p the total laboratory momentum of the dilepton. 
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For p > 5 GeV/c this error is of order 1%. Moreover, the directly produced J/ip 
[Y(1S)] and ijj' [T(2/3S f )] are expected to have the same production mechanisms 
and, therefore, very similar polarizations. As a consequence, the polarization ofj/ip 
[T(15)] from ip f [T(2/35)] can be considered to be almost equal to the polarization 
of directly produced J/ip [T(15)], so that, at least in first approximation, the two 
contributions can be treated as one. 

On the contrary, the J/i/j [T(15)] mesons resulting from XcJ [Xbj] radiative de- 
cays can have very different polarizations with respect to the directly produced ones. 
Directly produced P and S states can originate from different partonic and long- 
distance processes, given their different angular momentum and parity properties. 
Moreover, the emission of the spin-1 and always transversely polarized photon nec- 
essarily changes the angular momentum projection of the qq system in the P — >> S 
radiative transition. As a result, the relation between the "spin-alignment" of the 
directly produced P or S state and the shape of the observed dilepton angular distri- 
bution is totally different in the two cases: for example, if directly produced J/^, Xci 
and Xc2 all had "longitudinal" polarization (angular momentum projection J z = 
along a given quantization axis), the shape of the dilepton distribution would be of 
the kind 1 — cos 2 $ for the direct J/^, 1 + cos 2 ^ for the J/ip from Xci and 1 — § cos 2 $ 
for the J/ip from Xc2- While for directly produced S states — 1 < A# < +1, for those 
from decays of Pi and P^ states the lower bound is —1/3 and —3/5, respectively. 
More detailed constraints on the three anisotropy parameters A#, A^ and A#^ in the 
cases of directly produced S state and S states from decays of Pi and P2 states can 
be found in Ref. [26J Figure 3 of that work shows that the allowed parameter space 
of the decay anisotropy parameters for the directly produced J/ip [T(1S)] strictly 
includes the one of the S-states from P2 decays, which, in turn, strictly includes the 
one of the ^-states from Pi decays. 

The feed-down fractions are not well-known experimentally. In the charmonium 
case, the Xc-to-J/^ and Xc2-to-x c i yield ratios have been measured by CDFEZEH 
in the rapidity interval \y\ < 0.6, with insufficient precision to indicate or exclude 
important pt dependencies. The px-averaged results, 

R(Xa) + R(Xc2) =0.30 ±0.06, 

(6) 

R(xc2)/R(xd) = 0.40 ±0.02, 

where R(xci) an d R{Xc2) are the fractions of prompt J/i/j yield due to the radiative 
decays of Xci and Xc2, effectively correspond to a phase-space region (low pt and 
central rapidity), much smaller than the one covered by the LHC experiments. 

CDF also measured^ the fractions of Y(15) mesons coming from radiative 
decays of IP and 2P states as, respectively, R(xbi) + R{Xb2) — (27 ±8)% and 
R(Xbi) + R(x'b2) ~ (11 ±5)%, for > 8 GeV/c and without discrimination between 
the J = 1 and J = 2 states. These results tend to indicate that the contribution of 
the feed-down from P states to T(1S) production is at least as large as in the corre- 
sponding charmonium case, even if the experimental error is quite large. The same 
indication is provided with higher significance by the T polarization measurement 
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of E866 14 , at low px, as discussed below. 

Using available experimental and theoretical information, we can derive two illus- 
trative scenarios for the polarizations of the charmonium and bottomonium families. 
Figure IT] illustrates how the CDF measurement of prompt- J ftp polarizatiorP can be 
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Fig. 7. Direct- J/ip polarizations (A#) extrapolated from the CDF measurement of prompt- J/ip 
polarization (in the helicity frame), using several scenarios for the Xc polarizations. 



translated in a range of possible values of the direct- J/ip polarization, using the 
available information about the feed-down fractions and all possible combinations 
of hypotheses of pure polarization states for Xci and Xc2- The feed-down fraction 
is set to 0.42, two standard deviations higher than the central CDF value (Eq. [6|; 
using 0.30 simply decreases the spread between the curves. The R{Xc2) / R{Xci) ra- 
tio is set to 0.40; changes remaining compatible with the CDF measurement give 
almost identical curves. In the scenario in which Xci and Xc2 are produced with, 
respectively, J z = and J z = ±2 polarizations the CDF measurement is seen to be 
described by partial next-to-next-to- leading order (NNLO*) CSM predictions for 
directly produced ^-statefMHl The validity of this J ftp polarization scenario can 
be probed by experiments able to discriminate if the J/ip is produced together with 
a photon such that the two are compatible with being Xci or Xc2 decay products. 
Such dilepton events, resulting from Xc decays, should show a full transverse polar- 
ization (A^ cl = A^ c2 = +1), while the directly produced J/ip mesons should have a 
strong longitudinal polarization (A^ ir ~ —0.6). 

We will base our second scenario, for the bottomonium family, on the precise 
and detailed measurement of E866 14 , shown in Fig. [8^l. This result offers several 
interesting cues. It is remarkable that the T(26') and T(3S) are found to be almost 
fully polarized, while the T(IS') is only weakly polarized. The most reasonable 
explanation of this fact is that the fraction of T(IS') mesons coming from Xb decays 
is large and its polarization is very different with respect to the polarization of 
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Fig. 8. The E866 measurement of Y polarizations in the CS frame as a function of p^ (a), the 
deduced ranges for the fraction of T(1S) mesons coming from Xb decays (b) and the deduced range 
of their possible polarizations (c). A systematic uncertainty of ±0.06 is not included in the error 
bars of the data points in (a). The error bars in the derived lower limit for R(xb) (b) reflect the 
uncertainty in the A# measurements, assuming that the global systematic uncertainty affects the 
Y(15) and T(2S) + T(3S) measurements in a fully correlated way. The lower limits for R(xb) an d 
Atf(Yfromxb) depend on the ratio R(xb2) I R(Xbi) ? f° r which three different values are assumed. 



the directly produced T(1S). In fact, in the assumption that all directly produced 
S states have the same polarization, we can translate the E866 measurement into 
a lower limit for the feed-down fraction R(xb) from P states, summing together 
lPi, IP27 2Pi, 2P2 contributions. We assume that the T(26') + T(36') result has a 
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negligible contamination from x'b ~~ ^ T (2*5)7 decays and, therefore, provides a good 
evaluation of the polarization of the directly produced S states (a conservative 
assumption for this specific calculation). The lower limit for R(xb) corresponds to 

the case J z (xbi) = Jzix'bi) = ±:L > J z(Xb2) = JziXw) = °> in wmcn tne T(l£) 
mesons from Xb decays have the largest negative value of A# . The result, depending 
slightly on the assumed ratio between P 2 and Pi feed-down contributions, is shown 
in Fig. ||) as a function of px- More than 50% of the T(15) are produced from 
P states for (p T ) ~ 0.5 GeV/c, and more than 30% for (p T ) ~ 2.3 GeV/c. These 
limits are appreciably higher than the value of the feed-down fraction of J/ip from 
Xc measured at similar energy, low px and mid rapidity^. We remind that we have 
obtained only a lower limit (no upper limit is implied by the data), corresponding 
to the case in which Xbi and Xb2 are always produced in the same very specific 
and pure angular momentum configurations. Any deviation from this extreme case 
would lead to higher values of the indirectly determined feed-down fraction. 

The E866 measurement data also set an upper limit on the combined polariza- 
tion of Xbi and x&2- Figure^ shows the derived range of possible polarizations of 
T(15) coming from Xb- The upper bound, corresponding to R(xb) — 1 5 coincides 
with the measured T(15) polarization. The lower bound, slightly depending on the 
relative contribution of Xbi and x&2, is not influenced by the E866 data and corre- 
sponds to the minimum (px dependent) value of R(xb) represented in Fig. [8]o. The 
second strong indication of the E866 data is, therefore, that at low px the T(1S) 
coming from Xb decays has a longitudinal component in the CS frame larger than 
~ 30% (A# < 0.1), being ~ 60% (A# ~ —0.5) the maximum amount of longitudinal 
polarization that the T(1S) produced in this way is allowed to have. 

In the light of these scenarios it is clear that measurements of the polarization 
of the x states will be extremely important for an unambiguous understanding of 
the J/ip and Y(1S) polarizations. 

It has been shown in Ref. [26| that the angular momentum compositions of the 
X states produced in high energy collisions can be derived from the dilepton decay 
distributions of the daughter or T mesons, with a reduced dependence on the 
details of reconstruction and simulation of the radiated photon. This method is 
based on a particular choice of the quantization axes. Different frame definitions 
are in principle suitable for Xc and Xb polarization studies in hadronic collisions. 
We generically denote by V the charmonium and bottomonium 3 Si states, J/ip and 
T, and by x the 3 Pj states, Xcj and Xbj, with j = 1,2. Notations for axes and angles 
for the description of the x ~~ ^ ^7 decay are defined in Fig. [9j where z is the x 
polarization axis (for example the HX or CS axes, defined in the x res t frame). 
The traditional choice of axes, adopted in calculations and measurements (cited 
in Ref. [26]) of the full decay angular distribution for Xc mesons produced at low 
laboratory momentum, is represented in Fig. [lO^a), where the V polarization axis, 
z' is the V direction in the x res t frame. With respect to this system of axes, taking 
the polar anisotropy parameter A# as an example, all measurements will find, for 
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Fig. 9. Definition of axes and decay angles for x ~^ Vl- 
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where h 2 and g 2 are the fractional contributions of the magnetic quadrupole tran- 
sitions (electric octupole transitions for the j = 2 case have been neglected). The 
dilepton distribution in the x' ,y',z' coordinate system is independent of the x po- 
larization state. This choice of axes is suitable for measuring the contribution of the 
higher-order multipoles, but it does not provide information on the polarization of 
the x an d? hence, on its production mechanism, when the photon distribution is 
integrated out. An alternative definition, proposed in Ref. [26j enables the determi- 




nation of the x polarization in high-momentum experiments without the need of 
measuring the full photon-dilepton kinematic correlations. This definition, shown in 
Fig.fToFb), "clones" the x polarization frame, defined in the x res ^ frame, into the V 
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rest frame, taking the x" , y" , z" axes to be parallel to the x, y, z axes. As explained 
in Ref. |26j the dilepton distribution in this frame contains as much information as 
the photon distribution regarding the x polarization state: the two distributions are 
even identical when higher-order multipoles are neglected. 

The definition of the x, y, z axes (and, therefore, of the x h y" , z" axes) uses the 
momenta of the colliding hadrons as seen in the x rest frame, so that it requires, in 
general, the knowledge of the photon momentum. However, for sufficiently high (to- 
tal) momentum of the dilepton, the x and V rest frames coincide and the x" y", z" 
axes can be approximately defined using only momenta seen in the V rest frame. For 
example, if the x polarization axis (z) is defined along the bisector of the beam mo- 
menta in the x res ^ frame (CS frame), the corresponding z" axis is approximated by 
the bisector of the beam momenta in the V rest frame. The relative error induced by 
this approximation on the polar anisotropy parameter is |AA#/A# | = O [(AM/p) 2 ] , 
where AM is the x~V mass difference and p is the total laboratory momentum of 
the dilepton. Therefore, for not-too-small momentum this frame definition coincides 
with the frame defined in the measurement of the polarization of inclusively pro- 
duced / T mesons (CS or HX, for example). In other words, the measurement of 
the dilepton distribution at sufficiently high laboratory momentum provides a direct 
determination of the x polarization along the chosen polarization axis. This deter- 
mination is cleaner than the one using the photon distribution in the x res t frame, 
because it is independent of the knowledge of the higher-order photon multipoles. 

We remark that the general equations describing the x angular decay distri- 
butions depend on the exact definition of the quantization axes of x and V. An 
incorrect match between equations and frame definitions created confusion in the 
past, probably leading to wrong measurements^. While reading the recent LHCb 
papeJ^of Xc production in pp collisions we wonder if we might be in the presence of 
a similar situation. To calculate the effect of the unknown Xd and Xc2 polarizations 
on the calculation of the Xc and J/ip reconstruction and selection efficiencies, the 
authors reweight the simulated events assuming different Xd and Xc2 polarization 
scenarios. The formulas used for the angular distributions are those listed in the 
HERA-B paper on Xc production^, where the choice of the polarization axes is of 
the type represented in Fig. [lO[a). However, the axis definitions used in the LHCb 



paper seem to correspond to the convention represented in Fig. 10 ^b), considered 
in the high-momentum limit (certainly an excellent approximation in the LHCb 
case). First they define, to describe the prompt- J/ip decays, the angle Oj/^ as "the 
angle between the directions of the /i + in the J/tp rest frame and [of] the J/tp in the 
laboratory frame": that is, the chosen J/tp polarization axis is the centre-of-mass 
helicity frame. What perplexes us is the subsequent definition of the angles of the 
Xc —> ^ + ^~7 decay chain: "The Xc —> system is described by Oj/^ 

and two further angles, 6 Xc and (/>, where 6 Xc is the angle between the directions 
of the J ftp in the Xc rest frame and [of] the Xc in the laboratory frame [i.e., the 
ang le 8 of Fig. [9]". This means that referred to the J/ip helicity axis in the 
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laboratory, is taken as polar angle of the dilepton decays, while, to be consistent 
with the formulas used, a new angle should be introduced and defined as the 
angle between the directions of the /i + in the J/ip rest frame and of the J ftp in the 
Xc r est frame. The dilepton azimuthal angle 0, now defined as "the angle between 
the plane formed (...) [by] the Xc and J/ip momentum vectors in the laboratory 
frame and the J/ip decay plane in the J/ip rest frame", should rather be defined 
as "the angle between the J/ip decay plane in the Xc rest frame [the two leptons 
are collinear in the J/tp rest frame] and the plane formed by the Xc direction in the 
laboratory frame [being the helicity axis the chosen quantization axis for the Xc] 
and the J/ip direction in the Xc rest frame". 

We emphasize that such an inconsistency between axis definitions and formulas 
results in a completely wrong description of the angular distributions. The sentence 
"The angular distributions are independent of the choice of polarisation axis" may 
indicate a crucial misunderstanding: it is true that the angular distributions do not 
depend in form on the choice of the Xc quantization axis, but they depend drastically, 
also in form, on the choice of the J/ip quantization axis. Let us consider, for example, 
the dilepton distribution, integrated over the photon distribution, in the Xdi case. 
With the choice of axes made in the LHCb paper, and using the correct formulas, 
the cases of Xdi having helicity ±2 or would be observed, respectively, as a fully 
transverse (A# = +1) or dominantly longitudinal (A# = —3/5) J/ip polarization. 
Instead, the inconsistent formulas predict (erroneously, given the mismatch with 
the definition of the axes) an almost isotropic decay distribution (A# = +1/13), 
independently of the Xc2 helicity, as previously discussed. 

6. Polarization as an indication of sequential suppression 

Hypotheses on the suppression of Xc and Xb production in nucleus-nucleus collisions 
play a crucial role in the interpretation of the J/ip and T(15) measurements from 

gpg34|35|36^ RHICF|38|39|4U] and LHC |41|42|43| m termg of evidence Q f qua rk-gluon 

plasma (QGP) formation. The observation of the Xc and Xb suppression patterns 
in Pb-Pb collisions at the LHC could confirm or falsify the "sequential quarkonium 
melting" scenaridMl-^l and, therefore, discriminate between the QGP interpretation 
and other options. However, a direct observation of the Xc and Xb signals in their 
radiative decays to J/ip and Y(IS') is practically impossible in heavy-ion collisions, 
given the very large number of background photons produced in such events. 

The E866 scenario suggests an alternative method to determine the relative yield 
of P and S states by performing only dilepton polarization measurements. This pos- 
sibility is particularly valuable in the perspective of quarkonium measurements in 
heavy-ion collisions. A change of the observed J ftp and T(1S) polarizations from 
proton-proton to central nucleus-nucleus collisions would directly reflect differences 
in the nuclear dissociation patterns of S and P statesP^ Figure [ll] illustrates the 
concept of the method. The left panel shows an hypothetical R(Xc) pattern in- 
spired from the sequential charmonium suppression scenario, in which the Xc yield 
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disappears rapidly beyond a critical value of the number of nucleons participat- 
ing in the interaction (AT part ). This effect would be reflected by a change in the 
observed prompt- J/^ polarization. As shown in the right panel, according to the 
scenario presented in Fig. [7] the polarization should become significantly more longi- 
tudinal (in the helicity frame) after the disappearance of the transversely polarized 
feed-down contribution due to Xc decays. We are assuming that the "base" polar- 
izations of the directly produced S and P states remain essentially unaffected by 
the nuclear medium and are, therefore, not distinguishable from those measurable 
in pp collisions. A test of the sequential suppression pattern can, therefore, be made 
by comparing the prompt- J/ip polarization measured in pp (or peripheral nucleus- 
nucleus) collisions with the one measured in central nucleus-nucleus collisions and 
checking that this latter tends to the polarization of the directly produced S states, 
also determined in pp collisions through ip' measurements. 




350 400 



part '^part 

Fig. 11. A hypothetical variation of R(xc) (normalized to the pp value) with the centrality of 
the Pb-Pb collision (left) and the consequent variation of the prompt- J/ip polarization A# (right), 
according to the charmonium polarization scenario discussed in the text. 

The same method can be applied to the measurement of Xb suppression using 
T(1S) polarization. According to the E866 scenario (Fig. [8]), in pp (and peripheral 
Pb-Pb) collisions the T(IS') should be only slightly polarized, reflecting the mixture 
of directly and indirectly produced states with opposite polarizations. In central Pb- 
Pb collisions the T(16') would acquire the fully transverse polarization characteristic 
of the directly produced S states, indicating the suppression of the P states. 

We have estimated that about 30k prompt-J/^ and 10k T(16') signal events, 
with both leptons having > 5 GeV/c and an assumed background fraction of 
40%, would lead to a significant indication of the nuclear disassociation of the x 
states according to the scenarios we have considered. 

7. Summary 

Several puzzles affect the existing measurements of quarkonium polarization. The 
experimental determination of the J/iJj and T polarizations must be improved. 
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Measurements and calculations of vector quarkonium polarization should pro- 
vide results for the full dilepton decay angular distribution (a three-parameter func- 
tion) and not only for the polar anisotropy parameter. Only in this way can the 
measurements and calculations represent unambiguous determinations of the aver- 
age angular momentum composition of the produced quarkonium state in terms of 
the three base eigenstates, with J z = +1,0, —1. 

Moreover, it is advisable to perform the experimental analyses in at least two 
different polarization frames. In fact, the self-evidence of certain signature polariza- 
tion cases (e.g. a full polarization with respect to a specific axis) can be spoiled by an 
unfortunate choice of the reference frame, which can lead to artificial ("extrinsic") 
dependencies of the results on the kinematics and on the experimental acceptance. 

The angular distribution can be characterized by a frame-independent quan- 
tity, A, calculable in terms of the polar and azimuthal anisotropy parameters. This 
frame-invariant observable can be used during the data analysis phase to perform 
self-consistency checks that can expose previously unaccounted biases, caused, for 
instance, by the detector limitations or by the event selection criteria. The variable A 
also provides relevant physical information: it characterizes the shape of the angular 
distribution, reflecting "intrinsic" spin- alignment properties of the decaying state, 
irrespectively of the specific geometrical framework chosen by the observer. Extrin- 
sic dependencies on kinematics and acceptances are cancelled exactly, enabling more 
robust comparisons with other experiments and with theory. 

Stripped-down analyses which only measure the polar anisotropy in a single ref- 
erence frame, as often done in past experiments, give more information about the 
frame selected by the analyst ("is the adopted quantization direction an optimal 
choice?") than about the physical properties of the produced quarkonium ("along 
which direction is the spin aligned, on average?"). For example, a natural longitu- 
dinal polarization will give any desired A# value, from —1 to +1, if observed from 
a suitably chosen reference frame. Lack of statistics is not a reason to "reduce the 
number of free parameters" if the resulting measurements become ambiguous. 

Besides improving methodology aspects, more detailed and elementary informa- 
tion will have to be provided, by measuring separately the polarizations of directly 
and indirectly produced states. 

The forthcoming measurements of quarkonium polarization in proton-proton 
collisions at the LHC have the potential of providing a very important step forward 
in our understanding of quarkonium production, if the experiments adopt a more 
robust analysis framework, incorporating the ideas presented here. 

Quarkonium polarization can also be used as a new probe for the formation 
of a deconfined medium. This method, based on the study of dilepton kinematics 
alone, provides a feasible and clean alternative to the direct measurement of the 
X yields through reconstruction of radiative decays. With sizeable and T(15') 
event samples to be collected in nucleus-nucleus collisions, the LHC experiments 
have the potential to provide a clear insight into the role of the \ states in the 
dissociation of quarkonia, a crucial step forward in establishing the validity of the 
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sequential melting mechanism. 
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